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^ ABSTRACT 

i-Q ' We develop a systematic chiral perturbation expansion for the calculation 

of meson-exchange currents in nuclei and apply the formalism to nuclear axial 

/\ ' currents. We summarize the principal results of such a calculation to one loop 

C^ ' order on nuclear axial-charge transitions which provides a strong support to the 

conjecture of "chiral filter phenomenon" in nuclear medium. The use of heavy 
baryon chiral perturbation theory enables us to obtain a remarkably simple 
result valid next to the leading order in chiral counting. The dominant role of a 
soft-pion exchange in axial-charge transitions in heavy nuclei is confirmed. An 
important, alheit indirect, consequence of our result on the empirically observed 
enhancement in axial-charge transitions in heavy nuclei is pointed out. 



While meson-exchange currents in nuclei are fairly well understood in low-energy and 
low-momentum regime both experimentally and theoretically [l| , there remains the issue of 
understanding them from the point of view of QCD. This is not an academic issue since the 
purported high-energy electron machines in construction or in project are to probe nuclei 
for signals of direct QCD degrees of freedom in the deviation from exchange-current effects 
and for this, a contact with QCD at lower energy and highly nonperturbative regime will be 
clearly needed. At present the only applicable method to address this issue in nonpertur- 
bative regime is chiral perturbation theory (ChPT) ^. Chiral perturbation theory has had 
much success in Goldstone boson (vr, K) interactions]^, making nonperturbative QCD ac- 
cessible to laboratory phenomena. However incorporating baryons in the scheme has proven 
to be quite difficult, the main reason being that the standard power counting used in ChPT 
does not apply when baryons are involved. Even for the elementary vr-N interaction, the 
calculation becomes horrendously complicated, requiring approximations that are hard to 
control. Some progress has recently been made in describing by chiral perturbation expan- 
sion such processes as tt-N scattering [^, threshold pion production from nucleon [^ and 
nucleon's electromagnetic polarizabilities [^. Nonetheless a systematic chiral perturbation 
calculation of many-body nuclear properties remains a formidable, if not hopeless, task at 
the present stage of development. This technical difficulty explains in part the paucity up 
to date of work of this nature in the literature. Recently, however, this task was greatly 
facilitated by the formulation of heavy-baryon chiral perturbation theory (which we shall 
call HBChPT in short) by Weinberg [^j, Jenkins and Manohar Q and others g] along the 
line developed for heavy-quark effective theory B. 

The purpose of this letter is to summarize the principal results of a complete chiral 
perturbation calculation of the meson-exchange axial-charge operator to one-loop order that 
establishes rigorously the "chiral filter phenomenon" in nuclei conjectured a long time ago 
1 10 1 and given a partial justification recently by one of the authors [O] which states that 



whenever kinematically unsuppressed, soft-pion exchanges should dominate in electroweak 



processes in nuclei. As discussed in [10|, the chiral filter phenomenon occurs in nuclear axial- 
charge transitions and in nuclear Ml transitions (such as threshold radiative np capture). 
In this paper, we will focus on the former. The latter will be discussed in a separate paper 
in preparation |12|. Towner jl^ and Riska et al |14| have recently addressed a similar issue 
within, however, a phenomenological framework. We shall make comparison with their 
results at the end of the paper. 

As in [^, |l^] , we shall take the most general chiral Lagrangian consisting of nucleons 
and pions with all other degrees of freedom integrated out. (One can alternatively take a 
Lagrangian that also contains vector mesons and A's. We have satisfied ourselves that these 
elements do not alter our result significantly. This matter will be discussed in detail in a 
later publication.) The key observation made by the authors in [^, ^ is that the standard 
derivative expansion successful in low energy pion dynamics breaks down when baryons are 
involved, for the reason that time derivatives on the baryon field, typically of the order 



ruB ~ 1 GeV, are not small on the chiral scale set by A,^ ~ 1 GeV. The heavy-baryon 
formalism HBChPT circumvents this difficulty. To set up the formalism, one first redefines 
the velocity-dependent baryon field [^ -B„(x) = e^"^^^''"'"'^N{x) constrained to ^B^ = By 
where v^ is the four- velocity with f ^ = 1 and niB the mass of the baryon. This redefinition 
of the baryon field renders the relevant momentum scale of the baryon involved to be of the 
same order as the pion momentum. Consequently, HBChPT allows a chiral expansion in 
power of ((9/A^) and {d/niB) corresponding to the pion momentum scaled by A^ and the 
residual baryon momentum (when niBV^ is removed) scaled by m^, respectively. Expressed 
in terms of the B field, the chiral Lagrangian of pions and nucleons takes the formQ 
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where 2?^ = 9/^ + i^T ■ vf x DuTt and Dn = (1 -|- tt^/F^) ^d^ = D ^d^ are the "covariant 

TV 

derivatives" and S^ = i^5[^,"f'^] is the spin operator satisfying v^Sf,} = 0. The Ta in 
the quartic fermion term of Eq.(|l|) stands for all possible terms containing no derivatives 
allowed by Lorentz invariance and chiral symmetry and the ellipsis for higher-derivative and 
quark mass terms allowed by symmetries. In our case, higher derivatives acting on pion 
fields and quark mass terms play no essential role, so will be ignored. We will later make 
a specific reference to higher derivative terms involving baryons in connection with our 
renormalization procedure. The merit of this Lagrangian is that it provides a consistent 
chiral expansion for both mesons and baryons in terms of derivatives and quark masses 
P, R].[| Baryon momentum-dependent terms appear as higher order interaction terms in the 
chiral counting. The pair terms are suppressed at the leading order, appearing as "l/?n," 
corrections subsumed in the ellipsis in (|l[). 

Our task is to compute higher-order corrections in the chiral expansion parameter Q - 
where Q is the momentum or energy scale parameter probed by an external field, assumed 
to be small compared with the chiral scale A^ - to the leading soft-pion amplitude, de- 
noted Msofti that contributes in two-body effective currents responding to a slowly varying 
electroweak field J^, 

M=Msoft{l + 5 + 0{Q'')), 5r^0{Q^), n > 3. (2) 

There are no corrections of 0{Q) to 5 for the same reason that they are absent in nucleon- 
nucleon potentials Q. This chiral counting was established in jl^ for both electromagnetic 
and axial exchange currents. In what follows we will focus on the axial-charge operator 

* We can equally well rewrite this form in terms of the Sugawara field U = e'^^^'^'" frequently used in the 
literature |[|l. 

^The HBChPT as formulated is essentially equivalent to making static approximation on baryons fields, 
with the baryon velocity effectively conserved and with typical off-shell momentum of the baryon counted 
as of the same order as the pion momentum. 



as it is currently receiving considerable phenomenological attention [15|. We will use the 



axial current A^j^ obtained from the Lagrangian (|l|) by Noether construction, i.e., A^ = 

As emphasized recently by Weinberg p, |^, one should restrict the use of chiral 
perturbation strictly to irreducible graphs that are free of small energy denominators re- 
sponsible for binding. Exchange currents belong to this class of graphs. Now corrections 
to the soft-pion result can come from two sources: from contact four-fermion interactions 
involving derivatives and from one-loop graphs. We shall argue below that there should be 
no contribution from the former to the order we are considering. As for the latter, there are 
three classes of one-loop graphs to be calculated: The first (denoted as A) consists of one- 
pion-exchange current with one-loop renormalization of either the internal pion line or an 
external nucleon line or the ttNN vertex or the Jfj^nNN vertex; the second (denoted as B) 
consists of one-loop two-pion-exchange current with or without nucleon intermediate states; 
and the third (denoted as C) involves the four-fermion contact interaction F in (^. Most of 
the diagrams in A can be absorbed into the standard perturbative renormalizations of the 
pion mass, the nucleon mass and wavefunction and of the constants F^^ and qa- Since to 
one loop in the HBChPT there are no momentum-dependent corrections to these quantities, 
they are trivial and do not directly figure in our result. The only nontrivial contribution 
in this class comes from the J^ttNN vertex, the relevant diagrams of which are given in 
Fig.A(a-n). The class-B graphs are given in Fig.B(a-/i) and the class-C in Fig.C(a-6). 

At first sight, the number of graphs may look daunting but in the HBChPT with the 
Lagrangian (|l|) , many of these graphs do not contribute to the chiral order we are interested 
in. Some graphs vanish identically due to isospin symmetry (viz, Fig. Be). Other graphs 
such as Figs.Ca, Ai, Aj, Ak, Al, Am and An vanish because they are proportional to v ■ Sy 
which is zero. The graphs Cb, Ag, Ah, B/, Jig and B/i are proportional to S^, hence do 
not contribute to the time component since S^ ~ 0{Q/mB)- Thus we are left with only 
four graphs B(a, b, c, d) in the classes B and C and six graphs A(a-/) in the class A. We 
should note that while contact four-fermion interactions are very important in nuclear forces 
1^, P , they are suppressed not only at the leading chiral order as shown in [0] but also at 
higher orders. This feature is manifest in Weinberg's form of nonlinear chiral Lagrangian, 
the consequence of which constitutes one of the key elements in this work. This is in stark 
contrast with nuclear forces where numerous counter terms that are not readily available 



from experiments make systematic loop corrections highly problematic |16]. 

We shall take the momentum carried by the current to be zero and evaluate the 
two-body amplitudes corresponding to the exchange axial-charge operator. The HBChPT 
renders the loop integrals doable analytically. We note that even after the usual (wave- 
function, mass and coupling constant) renormalization, there are additional (logarithmic) 
divergences left over (i.e., L defined below). This is expected as we are dealing with a 
nonrenormalizable theory. This however poses no difficulty in our case. As is the stan- 
dard practice ||2|, these divergences can be absorbed in the coefficients of the counter terms 



that are next order in chiral expansion. The counter terms that absorb all the one-loop 
divergences in our calculation - that are subsumed in the ellipsis in (||) - are (in Weinberg 
representation) of the following form involving higher derivatives: 
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with ^2 = K2 + 24^(1 + 5g\)r], d[' = k[> + j^iSg\ - 2)r] and 4 

r] = j^ + r'(l) + ln(47r) — ln(r7T,^), r'(l) ~ —0.577215 and k's are finite constant counter 

terms.[] Now performing the calculation with dimensional regularization, we find 
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where q^ the four-momentum transferred from nucleon '1' to nucleon '2' carried by pions, 



KoiQ' 



-2 + ay, KiiQ^) = l-^y and i^2(Q' 



4 I (T^ I (7(3-(T-') 

9 "•" 6 "•" 12 



y with 2/ = In 



£+1 
cr-l' 



/ Ath ^ -I- O 

cj = J — Qj . Note that all the constants appearing here are physical ones which should be 
identified with experimental values. We have replaced at appropriate places the momentum 
q2 by the expansion scale parameter —Q^. The ratio of the spin-isospin matrix elements 
£ = ) ^M^ J'^lf-^^^-^l I figuring in the two-pion contribution is introduced for convenience. 

^ {iTiXT2q-{cri+a2)} & & J:' 

This ratio can be estimated in various nuclear models: In simple jj shell model, in Wigner 
supermultiplet model as well as in Fermi gas model of the nucleus, the ratio comes out to 
be l| so we shall set £ = 1 in the numerical estimates made below. As defined in (^), the 
quantity 6{1tt) comes from J^j^TrNN vertex renormalization for which only the six graphs 
Figs. Aa-Af survive to the chiral order we are calculating and 6(2-11) from the surviving 
two-pion exchange graphs. 

A close inspection shows that the terms involving d^ in (|3[) cannot arise from single 
vector-meson exchanges or other excitations that are lower than the chiral scale A^ ~ 1 

* rj and d's are singular and contain logarithms of mass. Here we briefly sketch our renormalization 
prescription. In calculating loop graphs, we use dimensional regularization. We encounter singular quantities 
in the form of 77 as given above. To remove them, we write a counter-term Lagrangian which is formally of 
the same chiral order as the one-loop graphs we compute. We adjust the coefficients of this counter-term so 
as to obtain a regular expression. The constants di so introduced contain two parts: one is proportional to 77 
and removes divergences and the other, the finite constants Hi which are to be determined (in principle) from 
experiments. The renormalization will be done at on-shell point for the nucleon and at zero-momentum for 
the pion. To one loop, the pion mass m^ and the decay constant _FV are independent of the renormalization 
point and can be taken from experiments. 

^ We would like to thank K. Kubodera for an invaluable help on this ratio. 
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GeV. Furthermore a constant term proportional to k\ (plus other counter terms that 
are in principle present even though no regularization is required) implies a zero-range 
interaction depicting the exchange of very massive degrees of freedom. Since we are to apply 
chiral perturbation expansion to only irreducible graphs while all reducible graphs are to be 
taken into account in calculating nuclear wave functions from a Schrodinger equation (or a 
relativistic generalization thereof) with a potential defined with the irreducible graphs and 
consequently the nuclear wave functions so obtained must contain short-range correlations, 
the consistency with the scheme requires that when embedded in nuclear matter, such a 
contact term be suppressed by nucleon-nucleon correlations. This invites us to drop the 
constant terms (or 6 function terms in coordinate space) Q. The remaining constant K2, 
which figures in the calculation of Ai^^ttNN vertex, can be fixed by the isovector charge 
radius of the nucleon, 

K2 = -];F^{r'^)X ~ -0.0856. (5) 

D 

This results because the A^ttNN vertex is related to the isovector Dirac form factor of the 
nucleon, F^ (t). This can be understood also by current algebra or in terms of vector-meson 
exchange ]T^]. 

One can have a rough idea of how large the chiral loop corrections can be by taking, 
in eq.(Q), Q ~ m^ as befits the scale involved in the chiral expansion. It comes out to be 

\6{Q^m^)\< 0.05. (6) 

To make a more quantitative estimate in nuclear matter, we have to go to coordinate space. 
It is in this space that short-range correlations mentioned above are most straightforwardly 
taken into account||J] Let us write the two-body axial-charge operator as 



M = {1 + 6soft)Msoft + Mioop (7) 
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0.0455. (8) 



In (0), we have separated out the long-range {0{Q'^)) contribution, denoted Sgoft, from 
one-loop corrections to the one-soft-pion exchange. The remainder is shorter-ranged and 

^Purists might object to this procedure by arguing that one has to calculate both nuclear forces and 
current matrix elements to the same order of chiral expansion. This we believe is not the right way of using 
chiral perturbation theory in nuclei. In fact, it is not a fruitful way of doing physics as a little thought would 
reveal that such a so-called "consistent" approach is doomed to fail. This "failure" should however not be 
construed as a failure of ChPT in nuclear physics as some people seem to argue. 

"in momentum space, this procedure roughly corresponds to subtracting constant terms from the ex- 
pressions of (kf). This is analogous to the constant subtraction one does to incorporate the Lorentz-Lorenz 
effect in vr-nucleus scattering. However in the present case, because of non-analytic terms, such simple 
prescriptions are not reliable. We have no choice but to go to coordinate space. 



hence combined with two-pion-exchange contribution into Mioop, representing the nontrivial 
part of one-loop corrections. The expression for Mioop is rather involved and so will not 
be written down explicitly in this paper. See |12] for details. In Fig. 2 is plotted the 



quantity 47rr^A^(r) where r = |ri — r2|. The axial-charge matrix element in nuclear medium 
calculated in fermi-gas model with a hard-core cut-off of d = 0.4 — 0.7 fm is given in Fig. 
3 as a function of matter density. Two key features to note in these results are: First, the 
loop corrections are very small for r > 0.6 fm compared with the soft-pion term consistent 
with the chiral filter argument and second the density dependence of the loop correction is 
weak. More precisely, for the reasonable hard-core cut-off of d = 0.5 fm, the ratio of loop 
correction over soft-pion R ~ 0.067 for p = 0.5po and ~ 0.089 for p = pQ. Even at nuclear 
matter density, the loop correction represents less than 10% of the soft pion result. 

To summarize, we have shown that the loop corrections to the soft-pion exchange 
axial-charge operator can be easily calculated in the HBChPT that provides a consistent 
chiral expansion. In this formalism, nucleon-antinucleon pairs are suppressed to next to 
the leading order. They can only contribute at higher chiral orders. Here we focused on 
the axial-charge operator but the same calculation can be done with no greater difficulty 
for the space component of the vector current, e.g., the Ml operator relevant for the elec- 



trodisintegration of the deuteron where the soft-pion effect is even more spectacular [12|. 
An important outcome of our calculation to one loop is that the chiral filter mechanism 
is robust, the soft-pion term playing a predominant role with one-loop chiral corrections 
remaining in the noise at small momentum transfers. It is hopeful that it will survive higher 
loop corrections. It is pleasing that our result eq.(|^) is totally free of unknown parame- 
ters inherited from more massive degrees of freedom once it is accepted that short-range 
correlations are operative in the transition matrix elements. 

We should point out one important, though indirect, consequence of our result. Since 
the suppression of higher-order chiral corrections established in this paper is likely to persist 
independently of the environment into which the operator is embedded, it is natural to 
conclude that higher-order chiral corrections cannot generate the apparently significant 
density dependence of the meson-current enhancement in heavier nuclei observed in nature 
|jl5[. The present calculation strongly suggests that the origin of the apparent density 
dependence of the axial charge transitions observed in Nature lies outside of higher-order 
chiral corrections, thus pointing to the possibility that the basic mechanism for the enhanced 



axial charge in heavy nuclei is indeed, as proposed in |17|, the scaling of hadron masses and 
quark condensate in dense medium, which is an intrinsic vacuum property |18]. Since the 
argument is based simultaneously on chiral symmetry and scale anomaly of QCD, it is 
perfectly consistent with the chiral expansion. We note that an empirical support for the 
scaling notion comes from a recent experiment on ^^B{p,j^) at 200 MeV [ff^]. 

This conclusion raises the question as to how the phenomenological Lagrangian ap- 



proach of 1 13, 14 1 can be understood in terms of the result we have obtained here. In par- 



ticular, since the HBChPT relies on the suppression of pairs, it is natural to ask whether 

6 



or if so, how the pair terms involving heavy mesons of [13, |lj] can be interpreted. The 
answer to this question that we propose is that at least the main pair term that figures in 
1 13, 14 1, namely the a exchange, corresponds to replacing the mass of the nucleon niB by the 



medium quantity m^ in the sense suggested in |18], with the vector meson mediated pair 
terms suitably suppressed by short-range correlations. A similar suggestion has been made 



in [pOl. This means that while two-body pair term in |13, 14 1 renormalizes the mass in the 



single-particle operator, three-body operators involving pairs will be needed to renormalize 



the two-body operator in the way it figures in |1?]. It would be interesting to verify this by 



an explicit calculation with phenomenological Lagrangians. 
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FIGURE CAPTIONS 

Fig. 1 

• A: The class-A one-loop graphs that renormalize the J^nNN vertex con- 
tributing to one-pion exchange axial-charge operator. The cross represents 
the axial current, the solid line the nucleon and the dotted line the pion. 
Only the graphs (a), (b), (c), (d), (e) and (/) survive to contribute. 

• B: The class-B one-loop graphs (a)-(/i) for two-pion-exchange axial-charge 
operator. Only the graphs a, b, c and d survive. 

• C: The class-C one-loop graphs (a)-(6) involving four-fermion interaction. 
Both graphs do not contribute to the chiral order considered. 

Fig. 2 

• Two-body axial-charge operator as a function of the separation distance 
obtained by putting -fi x T2 (cJi -|- CJ2) • r = {ti + T2) ai x (T2 • r = 1. Solid 
line is 47rr^(l -|- 6soft)M.soft and dotted line is ATir"^ Mioop- 

Fig. 3 

• The ratio of the matrix elements R = {Aiioop)/{i^ + Ssoft)Msoft) as a 
function of p/po for various hard-core cut-off d in the cut-off function 6(r — 
d). 



